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Abstract 

This paper is the second part of a series of two papers dealing with bulking: 
a way to define quasi-order on cellular automata by comparing space-time di- 
agrams up to rescaling. In the present paper, we introduce three notions of 
simulation between cellular automata and study the quasi-order structures in- 
duced by these simulation relations on the whole set of cellular automata. Var- 
ious aspects of these quasi-orders are considered (induced equivalence relations, 
maximum elements, induced orders, etc) providing several formal tools allowing 
to classify cellular automata. 
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1. Introduction 

In the first paper Q , we have developped a general theory of bulking aimed 
at defining quasi-orders on cellular automata based on the idea of space-time 
rescaling. The present paper focuses on three instances of such quasi-orders 
and uses them as classification tools over the set of one-dimensional cellular 
automata. 

Classifying does not make sense without additional assumptions (some cri- 
teria of classification). If in Wolframs papers 39] these criteria were implicit 
and informal, several classifications with explicit and formal criteria have been 
since proposed 0, ■ Usually, the criteria are those of dynamical systems 
and consist in a finite list of qualitative behaviors. Our approach here is dif- 
ferent: we don't define any a priori list of behaviors. Instead, we consider a 
simulation relation (a quasi-order) which tells when some cellular automaton 
is able to reproduce the behavior of another. The criterion of classification is 
then the definition of the quasi-order. Our central thesis is that, when it comes 
to apprehend the great variety of behaviors in cellular automata, the language 
of orders (equivalence classes, chains, ideals, maximal elements, distance to the 
bottom, etc) is more adapted than a finite list of monadic predicates (of the 
form "having property P" for some P). 

In this paper, we introduce three quasi-orders. They are all defined accord- 
ing to the same scheme developped in the companion paper Q: some local 
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comparison relation up to spatio-temporal rescaling. They only differ in the 
local comparison they use, which are based on the two following basic notions: 

• the injection of a small system (A) into a larger one (£?), 

• the projection of a large system (B) onto a smaller one (A). 
The three quasi-orders can be defined informally as follows: 

• an injective simulation of A by B, denoted by A =4i B, is an injection of 
some rescaling of A into some rescaling of B; 

• a surjective simulation of A by B, denoted by A =^ s B, is the projection 
of some rescaling of B onto some rescaling of A; 

• a mixed simulation of A by B, denoted by A =^ m B, is the injection into 
some rescaling of B of some C that projects onto some rescaling of A. 

In the context of cellular automata, the two notions of local comparison 
above (injection and projection) translate into the following. The first notion is 
the subautomaton relation {B obtained from A by forgetting some states) and 
the second one is the quotient relation (B obtained from A by identifying some 
sates) . 

The subautomaton relation was already introduced in [8| and its importance 
in cellular automata is illustrated by theorems 4 and 5 of that paper. 

The quotient relation can be seen as a particular case of the notion of factor 
in dynamical systems theory and symbolic dynamics (homomorphism between 
shift-commuting continuous global maps, see |21|). More intuitively, the quo- 
tient relation is a means to extract coarse-grained information (A) from a com- 
plex system (B) (see (l7l|). For instance, the metaphor of particles moving in a 
stable background used in the literature of cellular automataJ3| follows this idea: 
some information (e.g. the phase of the background) is hidden by identification 
of states. However, our definition of quotient requires that both B (the original 
system) and A (the one obtained after identification of some states) are cellular 
automata. 

We study these orders with several points of view and aim at understanding 
their structure as well as showing that they suitably capture many classical 
properties or phenomena of cellular automata. 

For instance, concerning the phenomenon of universality, we show that orders 
=4i and =4 m have a maximum, that classes of CA having Turing-universality can 
be obtained by simulating (in a way closed to Smith III [33] ) a universal Turing 
machine, and that such a class is not necessarily at the top of the order. 

As another example, we show that many global properties of cellular au- 
tomata as dynamical systems (reversibility, sensitivity, expansivity, etc) or cel- 
lular automata as computational devices (ability to simulate a Turing head, or 
to propagate some signal) characterize an ideal or a filter in our orders. 

Overview of the paper. Section 11.11 introduces three different comparison rela- 
tions which are three different instances of the bulking theory developed in the 
companion paper Q . Section [2] sets the definitions of these three notions of 
simulations and establishes some of their basic properties. Section [3] studies the 
'bottom' of each of the three quasi-orders induced on CA, i.e. CA or classes of 
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CA of least complexity. Section 2] focuses on the order structure with respect 
to various classical properties of CA, and from a computability point of view. 
Then section [5] explores the set of CA at the 'top' of these quasi-orders: uni- 
versal CA. Once again, the point of view is both structural and computational. 
Finally, section [6] is devoted to the construction of noticeable induced orders 
(like infinite chains), and the study of how simple families of CA spread over 
these quasi-orders. 

1.1. Definitions 

In this paper, we adopt the setting of one-dimensional cellular automata 
with a canonical neighborhood (connected and centered). A cellular automaton 
(CA) is a triple A — (S, r, f) where: 

• S is the (finite) state set, 

• r is the neighborhood radius, 

• f : S 2r+l — > S is the local transition function. 

A coloring of the lattice Z with states from S (i.e. an element of S z ) is 
called a configuration. To A we associate a global function G acting on configu- 
rations by synchronous and uniform application of the local transition function. 
Formally, G : S z -> 5 Z is defined by: 

G(x) z = /(X z _ r , . . .,X z+r ) 

for all z G Z. Several CA can share the same global function although there 
are syntactically different (different radii and local functions). However we are 
mainly interested in global functions and will sometimes define CA through their 
global function without specifying particular syntactical representations. In 
addition, the Curtis-Heldund-Lyndon theorem [13| allows us to freely compose 
global CA functions to construct new CA without manipulating explicitly the 
underlying syntactical representation. 

When dealing with several CA simultaneously, we use index notation to 
denote their respective state sets, radii and local functions. For instance, to A 
we associate S_a, r_4 and f^. 

This paper will make an intensive use of PCS transforms defined in section 
4.2 of 8], but restricted to dimension 1. With this restriction, a PCS transform 
a has the form a — (m, r, T, s) where m and T are positive integers, s is a 
(possibly negative) integer and r is either 1 or — 1. 

For any CA A, we denote by A^ or more explicitly A^ m ' T ' T ' s ' the applica- 
tion of a to A, which is, according to notations of 8], a CA of state set and 
global rule: 

(ffl m , V T D m ) o a s o G T A o (ffl m , V T D,,,) -1 . 

To simplify notation we will use a shortcut for purely temporal transforms: for 
any CA A we denote by A 1 the CA 1/ 4< 1 > 1 > t >°>. Finally, as another special case, 
we denote by A^ the grouped instance of A of parameter n: it corresponds to 
the transform (n, 1, n, 0) (see [1] for a detailed exposition of grouping). 
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2. Canonical orders 

In this section we introduce the three bulking quasi-orders that are studied 
all along the paper. They are obtained by applying the bulking axiomatics 
developed in the companion paper Q to three 'canonical' relations between 
local rules of CA. 

Those thee 'canonical' relations are in turn based on two classical notions of 
morphism between local transition rules of CA: sub-automaton and quotient- 
automaton. As shown below, the three relations we consider are exactly the 
reflexive and transitive relations that can be defined by compositions of one or 
more such morphisms. 

2.1. From Three Local Relations to Three Bulking Quasi- Orders 

A sub- automaton is a restriction of a CA to a stable sub-alphabet. A quotient 
is a projection of a CA onto a smaller alphabet and compatible with the local 
transition rul^E Both define a kind of morphism between cellular automata: 

• A is a sub-automaton of B, denoted A'Q B, if there is an injective map 
l : Sj. — > Sb such that I o Ga = o I, where I : S\ — )• Sg denotes the 
uniform extension of l. We often write A C t B to make the map i explicit. 

• A is a quotient of B, denoted A < B, if there is a surjective (onto) map s 
from Sb to Sa such that s o G& — Ga ° s, where s : Sg — > S\ denotes the 
uniform extension of s. We also write A < s B to make the map s explicit. 

Relations C and < are quasi-orders (reflexive and transitive) and it is straight- 
forward to check that their induced equivalence relation is the relation of iso- 
morphism between cellular automata (equality up to state renaming) denoted 

by = 

It is also straightforward to check that C and < are incomparable (none of 
them is implied by the other one) . It is thus interesting to consider compositions 
of them. The composition of two relations R\ and R2 is the relation R\ ■ R2 
defined by 

Ri ■ R 2 = {(x, y) : 3z, (x, z) <E R\ and (z, y) g R 2 }. 

We denote by 1Z the set of relations obtained by (finite) composition of < and 
C. Any relation of 1Z is a priori interesting, but the following theorem justifies 
that we restrict to <, C and the composition < • C only. In the sequel < • C is 
denoted by <C and, as for C and <, we use the infix notation (A <C B). 

Theorem 2.1. 

1. any relation R g TZ is included in <H- (i.e., (A,B) g R implies A <C B) ; 

2. the transitive relations oflZ are exactly: <, C and <TZ . 



1 A quotient is a particular kind of factor, a classical notion in dynamical systems theory 
and symbolic dynamics 12011 
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Proof. We first prove that if A C • < B then A <C B, which is sufficient to 
prove assertion 1 by transitivity of C and of <. So consider A, B and C such 
that A C t C and C < s B. Then consider Q = s^ot^). We have G B (Q Z ) C Q z 
because 

s o G e (Q z ) = G c o s(Q z ) (because C < s B) 
= G C °T{S\) (by definition of Q) 
= T o G^(5^) (because .4 E t C) 

ci(s5). 

The CA ^ = (Q, rg, /g) is thus well-defined and by definition we have X \—B. 
Moreover, we have A <L-i os % because t _1 o s : Q — > A is well-defined and onto, 
and because 

o s o Gx = G^4 o o s 

since s o Gg = Gc ° s and o Gc = G.4 o i -1 over (t(A)) Z = s(Q z ) . Hence 
.A <SZ B and assertion 1 is proven. 

Given assertion 1 we have 1Z — {<, C, C • <, <C}. To prove assertion 2, it 
is thus sufficient to prove that C • < is not transitive. To do this, consider 
Sa = {0, . . . ,p — 1} with p prime, p > 5, and let a, ao, ai, 60, 01 be five distinct 
elements of S.4. Then consider .4, the CA with state set S4, radius f and local 
rule ] a defined by: 

!a\-i if x ^ a and y = ai, 

bi-i if x ^ a and y = b i: 

y + 1 mod p else. 

/^t depends only on two variables. Suppose now that there is some AC B with 
at least two states such that B <^ A. We will show that n must be one-to-one. 
Suppose for the sake of contradiction that there are distinct elements x and y 
in Sa such that ir(x) = n(y). Then, because Ja(*, a, z) = z + I mod p for any 
z, we have tt(x + 1 mod p) = 7r(y + 1 mod p) and more generally 

tt(x + i mod p) — ir(y + i mod p) 

for all i € N. So, supposing without loss of generality y > x, let k = y — x. We 
deduce from above that Tr(y) = n(y + jk mod p) for all j G N and, by elementary 
group theory, that 7r is constant equal to n(y) (because p is prime and k 7^ 0). 
This is in contradiction with the fact that tt has image Sb which has at least 
two elements. So n is one-to-one and B is isomorphic to A. Now consider C, 
the identity CA over state set Sq = {0, i}. Since C possesses 2 quiescent states 
and A has no quiescent state (straightforward from the definition of /a above), 
we have C0A. With the discussion above, we can conclude that CJ^-*~<A. 

However, we have B <E A because the states {ao, a\, bo, 61} induce a sub- 
automaton C of A which verifies B < s C where s : {ao, a\, bo, b{\ — > {0,f} is 
defined by s(di) = et s(bi) = 1. Assertion 2 follows since the relation <C is 
included in the composition of the relation C • < with itself. ■ 

Like C (already considered in Q), < and <C are quasi-orders on CA and 
therefore constitute natural candidates for the divide relation of bulking ax- 
iomatics (definition 8 of 0]). 
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Inspired by definition 14 of j8j, we now define 3 bulking quasi-orders using 
PCS transforms. 

Definition 2.1. B simulates A injectively, denoted A 4i B, if there exist two 
PCS transforms a and (3 such that A^ E B^ . 

We will occasionally use the notion of simulation by grouping introduced in 



24j and discussed in |8(: we denote by A B the fact that there are n and m 



such that A [n] E B [m] - This is a special case of the injective simulation above. 

Definition 2.2. B simulates A surjectively, denoted A 4s B, if there exist two 
PCS transforms a and (3 such that A {a) < B^ . 

Definition 2.3. B simulates A in a mixed way, denoted A 4m B, if there exist 
two PCS transforms a and j3 such that A^ <E B^' 

For each notion of simulation above, we say that the simulation is strong if 
the transformation a applied to the simulated CA is trivial: a = (1, 1, 1,0) so 
that A {a) = A. 

Theorem 2.2. (CA,4i), (CA,4 S ) and (CA,4 m ) are quasi-orders. 

Proof. We show that 4i and 4m correspond exactly to models of bulking 
developped in Q: the proof of theorem 15 of Q contains the case of injective 
simulation. The case of =^ m follows immediately (axiom (-B4) is straightforward 
and axiom (B5) is verified because <C contains C). For =^ s , the proof of each 
axiom is similar except for axiom (-B5). 

With or without axiom (-B5), theorem 10 of [Sj can be applied in each case 
and show the present theorem. ■ 

Lemma 2.1. Let < be any relation among C, < and <C. Then the following 
propositions are equivalent: 

• there exist two PCS transforms a and f3 such that A^ < B^' , 

• there exist a PCS transform a and an integer t such that A^ < B^\ 

• there exist a PCS transform (3 and an integer t such that A^ < B^ , 

Proof. We use the property of compatibility of relation < with respect to 
geometrical transforms (axiom B± of Q). The lemma follows from the following 
property: for any transform a, there exist a transform (3 and an integer t such 
that 

\/F : (F {a) ) ifi) = F [t] . 

If a = (m, t, z, 0), (3 can be chosen as the composition of (1, m, 0, 0), (1, 1, — z, 0) 
and (t, 1,0,0). ■ 

In the sequel, if ^ denotes a simulation quasi-order we denote by ~ the 
induced equivalence relation and by [A] the equivalence class of A with respect 
to ~. For instance, to 4i we associate the notations ~j and [A] i with the 
following meanings: 

A~ l B <^ A 4i B and B 4i A, 
[A], = {B:A~i B}. 
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We use similar notations for =<! s and =<! m . 

Before entering into details concerning various aspects of the three simulation 
relations defined above, we can already make a clear (yet informal) distinction 
between =^ and =^ TO on one hand, and =^ s on the other hand. For the two 
former, the simulation takes place on a subset of configurations and nothing 
can be said a priori about the behavior of the simulator outside this subset 
of configurations. For =4 S , however, the simulation occurs on any configuration 
and the simulator's behavior on any configuration is in some way affected by 
the simulation. Section l4~2l give several illustrations of this difference. 

2.2. First Properties 

We now establish a set of basic general facts about =4i, =4 S and =^ m while 
next sections of the paper focus on particular aspects. 

Theorem 2.3. Let A be any CA and =<! be any relation among =$j, =<; s and =4 m . 
Then it holds: 

1. there is some B G [A] having a quiescent state, 

2. there is some B € [A] with radius 1, 

3. _L =^ A where _L is the CA with a single state, 

4. A^ Ax B and A =4 B x A for any B. 

Proof. 

1. there exists some uniform configuration x and some t > 1 such that G Ax) = 
x. So A 1 has a quiescent state and it clearly belongs to [A]. 

2. A^' 1 ' 1 ' ^ admits a syntactical representation with radius 1 and clearly 
belongs to [A]. 

3. First, one always has _L < T A where it is the trivial surjection mapping 
each state of A to the single state of _!_. So assertion 3 is proven for =^ s 
and =4 m . Second, one has _L Cj B if B has a quiescent state where i is the 
trivial injection mapping the single state of _L to the quiescent state of B. 
Assertion 3 follows for =^ by assertion 1. 

4. We show only the first relation, the second being rigorously symmetric. 
First, one has always A Ax B where iri : Sa x Sg — > is the pro- 
jection over the first component. Second, if B has a quiescent state q, 
one has A Q,, A x B where t is the injection defined by l(x) = (x, q) for all 
x E S_a (the equality To Ga — (G ' A x Gg) ol is true over S^). If B has 
no quiescent state, just consider B l and apply the previous reasoning to 
obtain: 

A t QA* x B % = (A x B)' 

and thus A =^ A x B. 

m 

The three simulation quasi-orders are derived through bulking axiomatics 
from three different relations on local rules (see 12. ip . There is a priori no 
reason why the differences between local relations should extend to differences 
between the three simulation quasi-orders. The following theorem shows that 

=4 S and =^ m are nevertheless different and that =<;.; and =^ s are both strictly 
included in =<;,„. 
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Theorem 2.4. The relations and are incomparable (no inclusion in ei- 
ther direction). 

PROOF. We first show that there are CA A and B such that AC. B but A^B. 
Let A = a x (t _1 defined over states set S_a = {0, 1} X {0, 1} and let B be the 
CA of radius 1 defined over 5g = Sj, U {#} by: 



f B (x,y,z) = 



fA{x,y,z) tfx,y,zeS A , 
y else. 



One clearly has A C Id B. Now suppose A B. Without loss of generality 
we can assume that there are geometrical transforms a =< m, r, T, s > and 
P =< m', 1, V, > such that _4< Q > < w . But, by definition of B, there exists 
some state qo of B^ 3 ' (for instance # m ) which is left invariant by iteration 
of B^ whatever the context. Then n(qo) must be a state of A^ with the 
same property. This is impossible since either s ^ T or s =/= —T and thus some 
component of the future state of a cell of A^ is dependent of the state of a 
neighboring cell. 

Now we show that there are A and B such that A < B but A^lB and the 
theorem follows. Let A and B be the automata pictured on figure [T] A is a CA 
with two states, and 1, whose behavior is to reduce ranges of l's progressively 
until they reach size 1: at each time step the cells at each ends of a range of size 
3 or more are turned into state (only the right cell of range of size 2 is turned 
into 0). B has three states (0, 1 and 2) and has the following behavior: ranges 
of size 3 or more of non-zero states are reduced in a similar way by the two ends 
(states inside ranges are left unchanged), ranges of size 2 become an isolated 2 
(left cell becomes 2 and right cell 0), and ranges of size 1 become an isolated 
1. In a word, B reduces the size of non-zero ranges until size 1 but keeps the 
parity information at the end: an even range becomes eventually an isolated 2 
and an odd range becomes an isolated 1 (see figure [1]) . 

Formally, let 7r : {0, 1, 2} — > {0, 1} be the surjective function defined by 7r(0) = 
and 7r(x) = 1 if x ^ 0. Now let A be the CA of radius 2 and state set Sa = {0, 1, 2} 
with local rule: 



fA(x,y,z,t,u) 



1 if Tf(xyztu) — OHIO, 

2 if W(yztu) = 0110, 

z if Jf(yzt) — 111 and if Tr(xyztu) ^ OHIO, 

z if W(yzt) = 010, 

in any other case. 



Finally, let B be the CA with states set B = {0, 1}, radius 2 and local transition 
function 



f B (x,y,z,t,u) 



1 if yzt — 111 or yzt = 010 or yztu — 0110, 
else. 



By construction, we have A <,r B. Now suppose for the sake of contradiction 
that A =4i B and more precisely: 

A {m>yy.s') c B <m,M,0) 



8 




Figure 1: Behavior of A (left) and B (right). Time goes from bottom to top. 



where a =< ml ', r', t', s' > and j3 =< m, 1, t, > are suitable geometrical trans- 
forms. Let u = l m and u = m (it and v are particular states of A^) and 
consider J7 = 4>(u) and V = (U and V belong to Sg l ). The remaining of 
the proof below proceeds by a careful case analysis on U and V to obtain a final 
contradiction. The main technique is to consider specific orbits of A^ involv- 
ing u and v, and to derive constraints on their possible image by <j) involving U 
and V . 

Since configurations u and v are fixed points of A^ a \ so are U and V for 
B^'. Moreover, one can check from the definition above that the state is a 
'blocking state' for B: the half-configuration on the left of an occurrence of 
evolves independently of the half-configuration on its right. So, if U contains 
one or more zero's, then any configuration of B^ containing U 3 will contain at 
least one occurrence of U for ever (because it is the case for the configuration 
U): this is in contradiction with the fact that the orbit of a configuration of 
the form u vu A if does not contain any occurrence of u after sufficiently many 
iterations of A^ (because, whatever the value of to', u 4 represents in A an 
even-sized range of Is which is reduced until the last two Is are turned into a 
single 2 by case 2 of the definition of A). Hence we have W(U) = l m . 

From this we deduce that V = m because configurations of the form UJ VU n V i 
are transformed into configurations where a single cell is not in state V (just 
consider the orbit of u vuv u under .4^) and large ranges of non-zero states are 
always turned into large ranges of zero's under B. 

Finally, we have 7f(0(O m = m_1 l by considering the orbit of a config- 

uration of the form M vu 2 v^ under A^ and its counterpart of the form U VU 2 V U} 
under B^> (by the way, we also show that the shift parameter of transform a is 
0). Now, letting v! = m _1 f , we have on one hand the orbits of 2 configurations 
of the form u vu'u 2n v u and u um 2 "d u both leading to the same configuration of 
the form "vu'v^ under A^ , and on the other hand, the orbits of "V^u^U^V" 
and ^VXJ 2n V^ leading to different fixed points under B^ due to different parity 
of non-zero ranges: this is a contradiction since d> o A^ = B^ o 6. m 

3. Bottoms of the Orders 

This section focuses on the bottom of the orders. We have already seen 
(theorem !2.3p that T is a global minimum for the three quasi-orders considered 
here. In this section, we study CA that are at the lowest levels of the quasi- 
orders. Formally, the only CA at level is _L and a CA A is at level n + 1 for 
a quasi-order =<! if: 
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1. A is not at level n and, 

2. VB : B 4 A ^ B e [A] or B is at level i with i < n. 

The following theorem shows that some classical properties of CA correspond 
to classes at level 1. Recall that a cellular automaton is nilpotcnt if all initial 
configurations lead to the same configuration after a finite time. 

Theorem 3.1. Let =<! be a simulation relation among =<!j, =^ s and =4 m - Then 
the following CA are at level 1 (provided they have 2 or more states): 

1. the set of nilpotent CA, which is an equivalence class for ~, 

2. the set of CA which are periodic up to translation (A 1 o <j z = Id) which is 
exactly the equivalence class for ~ of the identity CA. 

Proof. 

1. Nilpotency is equivalent to the existence of a uniform configuration reached 
in a fixed finite time from any configuration. This property of phase space 
is clearly invariant by geometrical transforms and preserved by taking sub- 
automata or quotient automata. So any nilpotent CA is at level at most 
1. Moreover, the set of nilpotent CA forms an equivalence class. Indeed, 
for any nilpotent A, there is t such that A* is a constant function equal 
to some <7a. If we consider any nilpotent B with at least 2 states, there is 
m such that \Slg\ > \S^\ and t' such that £>* is a constant function equal 
to some 7]b. If we consider the geometrical transforms a —< 1, 1, t, > and 
(3 =< m, 1, i', >, then we have both A {a) C l B<^ and A {a) B^ if i is 
such that i(q a ) = qb and 7r is such that ir(x) = q a •<=>■ x = qb- 

2. Any CA which is periodic up to a translation is by definition equivalent 
to some identity CA and two identity CA with different state set are also 
clearly equivalent. Moreover, all such CA are at level 1 because the property 
of being periodic up to translations is preserved by geometrical transfor- 
mations and by taking sub-automata or quotient automata. ■ 

In the remaining part of this section, we will study two families of cellular 
automata with respect to the quasi-orders: a subset of additive CA and products 
of shifts. Our goal is to show that at (almost) each finite level there are infinitely 
many incomparable classes (theorem 13.31 and corollary 13.11 below) . 

3.1. Additive Cellular Automata 

The bottom of the quasi-order was studied in [24j]. The main result is 
the existence of an infinite familly of mutually incomparable CA at level 1: the 
familly of CA Z p with p a prime number and where Z p is a CA of radius 1 and 
state set {0, . . . ,p — 1} defined by the following local rule: 

<5z P ( x > V, z ) = x + V + z m °d P- 

There are strong connections between ^p and =<^ and in fact the set of CA 
at level 1 are the same for these two quasi-orders. 

Lemma 3.1. If A is at level 1 for ^p then A is at level 1 for =<!,. 
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Proof. If B =4i A then by lemma l2~Tl there is some integer t and some transform 
/? such that B^' C A™ . By theorem 12.31 we can suppose that A has radius 1 
so B'P' has radius 1. Since A is at level 1 for then either B^' € [A] a or 
B^ G [_L] n . We deduce that either B G [A^ or S £ [1^. Hence .4 is at level at 
most 1 for =4i and it cannot be at level since it is not in [_L] i = [J-]q. I 

The previous lemma is not enough to show that the CA {Z p j^ with p prime 
are mutually ^-incomparable because several equivalence classes for sCp can be 
included in a single class for =<!j. However we are going to show that this familly 
is a set of mutually incomparable CA for the three quasi-orders introduced 
abovfl Moreover, for ^rj and =4i, they are all at level 1. The proof relies on 
the following result already used for the case of 

A CA (S, r, f) is LR-permutative if the two following functions are bijections 
for all oil, . . . , ai r : 

• x H> /(oi, . . . , a2r, x) and 

f(x,ax, ■ ■ ■ ,a 2r ). 

Theorem 3.2 ([22;]). Let p be a prime number and t > 1. Then we have: 

1. Zp^ is LR-permutative; 

2. ifAQ Z p ^ then p divides \S_a\. 

To take into account use of < in simulation we will use the following lemma. 

Lemma 3.2. If B is LR-permutative and A <] B then \Sa\ divides \Sq\. 

Proof. To simplify notations, we suppose that B is of radius 1 (the proof 
works the same way for higher radii). Suppose A ^ n B. By surjectivity of n, it 
is sufficient to show that tt is balanced, i.e. such that for all x, y G Sa- 

\{e : 7r(e) = x}| = \{e : n(e) = y}\ . 

Consider any x, y € SU- Let a,b G Sjg be such that ir(a) = x and ir(b) = y 
and consider any c S By R-permutativity there is d G Sjg such that 

/e(a, c, d) — b. Now for any a' G Sb such that 7r(a') = 7r(a), we must have 
7r(/g(a', c, d)) = 7r(6) because .A < T £>. Moreover, by L-permutativity, a' H> /e(a', 
is one-to-one which proves: 

|{o : Tr(a) = 4| < |{6 : n(b) = y}\ . 

The balance of ir follows by symmetry. ■ 

The results above are the key ingredient of the following theorem. 

Theorem 3.3. Let =<! be any relation among =^j, =^ s and =^ m . Let p and q be 
two distinct prime numbers. Then we have: 

1- Zp ^ Zq 



2 The proof we give here was suggested by E. Jeandel. 
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2. Z p is at level 1 for 

Proof. 2 follows immediately from lemma [57X1 and the fact that Z p is at level 
1 for ^p (corollary 2 of |22[). To prove assertion 1 it is enough to prove 
Zp^frmZq- Suppose for the sake of contradiction that Z p =^ m Z q , or equiva- 
lent^ by lemma 12.11 that there are a CA A, a transform a and an integer t 
such that Zp^ <\ A.Q Z q ^. Then, combining lemma 15721 and theorem 13.21 we 
deduce that the number of states of Z p ^ is a power of q which contradicts the 
fact that p and q are two distinct primes. ■ 

3.2. Products of Shifts 

We will now study products of shifts in order to show that there are infinitely 
many incomparable CA at any finite level greater than 3 for any of the three 
simulation quasi-orders of the paper. 

We denote by a n<z the translation CA with n states {1, . . . , n} and transla- 
tion vector z defined by: 



We then consider cartesian products of such CA. Since cr„. z x a v z = a nPtZ , we 
can focus on considering cartesian products where all vectors are distinct. 

The next lemma shows that the structure of product of translations is pre- 
served when taking sub-automata and quotient-automata. 

Lemma 3.3. Let B = Y\i=i <7 n i ,z i (with Zi all distinct) and suppose A is such 
that A <n B then A is isomorphic to J\j=i °V»'. ,z» where 1 < ij < n and 



Proof. The lemma is straightforward if we replace <C by C. So it is enough 
to show that it is also true when replacing <C by <. Suppose that A < n B. 
The idea of the proof is to show that tt must be 'compatible' with the product 
structure: it forgets some components and keeps others but never introduces 
any kind of 'correlation' between them. So let i be such a component in B 
(1 < i < p). Consider any pair of states q, q' 6 S& such that qi = q[ (where qi 
or q\ denotes the projection on the ith component). Denote by g+ and q' + the 
states obtained from q and q' by changing their ith component in the same way 
((q+)i — (q' + )i). Since the Zi are distinct, one can build two configurations c 
and c' of B such that: 

• c(z) = c'(z) for all z ^ 0, 

• c(0) = q and c'(0) = q', 

• B{c) = q+ and 23(c')o = q' + . 

If 7r(q) = Tr(q') we have 7r(c) = 7f(c') so ir(q+) = Tr(q' + ). The same reasoning can 
be done starting from q + and q' + so we have: 



Hence, if there exist two states e and e' with the same image by tt and 
which agree on all components except component i, then values ei and e[ can 
be exchanged in the ith component of any state without affecting its image by 




2<n' h <n t] . 
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7r. In such a case, wc can consider the CA C obtained from B by identifying a 
and e- in the ith component. More precisely, C is of the form 

with n- = n, — 1 and = rij for any j 7^ i. Then we have .4 < C < S. Applying 
this reasoning iteratively, we finally have A <! 9 Co <! B where Co is of the form 
n,=i zi- where 1 < ij <n and 2 < n-. < rij, (component reduced to 1 state 
during one step of the process can be eliminated) and g is such that changing 
the value of any component of any state of Co will change its image by g. Now 
suppose for the sake of contradiction that g is not injective. Then there are 
states q and q' of Co such that g(q) = g(q'). Let c = q and d be equal to c 
except on position where it is in state q' . By hypothesis, at any position z, 
Gc (c) and Gc a (c') must be in states having the same image by g. But since Co 
is a product of translations with distinct vectors, there must be some position z 
where Gc ( c ) an d Gc ( c ') are m states which differ on one component only: this 
is in contradiction with the hypothesis on g. Hence g is injective and therefore 
A = C . M 

We will now study the effects of geometrical transformations on CA which 
are products of shifts. Of course, the translation vectors involved in such an 
automaton can be altered by geometrical transformations. If A is a product 
of translations with vectors z\ < . . . < z a , we denote by \ {A) the following 
characteristic sequence (provided a > 3): 

\ z 2 - z\ z 2 -zi) 

The purpose of the following theorem and lemma is to establish that the char- 
acteristic sequence gives a simple way to compare any pair of products of shifts 
in the three quasi-orders. 

Theorem 3.4. Let =<I be a relation among =4i, =^ s and =4 m . Let A be a prod- 
uct of a > 3 translations with distinct vectors and with characterisitic sequence 
X (A) = (ai, . . . , a a -2)- If B =4 A then B is equivalent to some C which is a 
product of a subset of b translations of A. Moreover, we have the following 
properties: 

1. ifb = a then C has the same characteristic sequence than A; 

2. ifb = a—l and b > 3 then the characteristic sequence of C has one of the 
following form: 

• (ai,...,ai_i,a i+ i,...,a _2) 

. fa 2-1) 

( "2-1 «a-2-l \ 

ai-1 J 

3. if the characteristic sequence of C is not \ (A) then A/{C. 
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Proof. Let z\ < z^ < . . . < z a be the ordered list of translation vectors of A. 
Since B ^ A, there is some C equivalent to B and some integer t > 1 such that 
C <C „4i'l (by lemma l2~Tj) . We deduce from lemma [5751 that C is isomorphic to 
a product of translations whose vectors are a subset of the familly (z^ since A 
and .Al'l have identical translation vectors, C must have the same characteris- 
tic sequence than A if it has the same number of translation vectors. When 
b = a — 1 and b > 3, it is straightforward to check that the three possible forms 
of the characteristic sequence of C correspond to the case where the missing 
vector is z,-, zi and z\ respectively. To prove the last assertion of the theorem, 
it is sufficient to check that for any transform a of the form < m, l,mt,mz > 
(we can restrict to such transforms by lemma IT27T1 C and C' a ' are products of 
translations with the same characteristic sequence because each vector Zi of C 
becomes tz\ + z in C' a '. ■ 

The next lemma gives canonical members of the equivalence classes of prod- 
ucts of shifts. 

Lemma 3.4. Let^ be the equivalence relation induced by any of the quasi-order 
=^ s and =4 m . Consider any t 0, any z and any product of translations of 
the form A — Yii=i "n,, tz t +z> ■ Then we have: 

A~ Y[ (72,z t 
l<i<p 

Proof. Let B = nf=i °"2,z 4 and let m = maxn^. It is straightforward to check 
that A < £<"M> mt > mz > and',4 C B^ m ^ mt ^ mz \ and also that B^-M.*) < A and 
B^'WcA ■ 

Theorem l3.4l and lemma [3~41 give a complete characterisation of the position 
of products of shifts in the quasi-orders considered in this paper. We will use it 
later in section 14.11 but we now state the main result of this section concerning 
levels at the bottom of the quasi-orders. 

Corollary 3.1. Let ^ be a relation among =4i, =4 S and =4 m . For any n > 3, 
there are infinitely many incomparable CA at level n for =<; . 

Proof. We have shown in theorem 13.11 that translations CA are at level 1. 
Lemma 13.31 together with lemma 13.41 show that a product of two translations 
(with distinct vectors) is at level 2. By theorem 13.41 we conclude that any 
product of n translations with distinct vectors is at level n and two such CA are 
incomparable if they have different characteristic sequences provided n > 3. I 

4. Structural properties 

In this section we study in various ways the order structures induced by the 
simulation relations defined above. 
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4-1- Cartesian Products and Lack of (Semi-) Lattice Structure 

The Cartesian Product of cellular automata is not a neutral operation from 
the point of view of the three quasi-orders of the paper. For instance, there are 
CA A, B which are equivalent but such that Ax A is not equivalent to B x A 
(it is sufficient to take two shifts with different translation vectors). 

The next theorem shows however that some simulations by Cartesian prod- 
ucts of CA can be transposed to components of the product. 

Theorem 4.1. Let A be a CA with two states and let =4 be a simulation relation 
among =4i, =^ s and =<; m . For any B and C, if B x C strongly ^-simulates A then 
either A ^ B or A 4 C. 

Notice that the 'two states' and 'strong simulation' hypotheses are both 
important and related. The theorem doesn't hold without such hypotheses: 
take two shifts with different vectors for B and C and choose A= B x C. 

Proof (of the theorem) . Let Sa = {ai, 0-2}- First, we show that A Q t B x C 
implies either A Q B or A C C which is sufficient to prove the theorem for =^ and 
=4 m . Since i(a\) ^ t(a 2 ) we have either ir\{i(a{)) ^ 71-1(1(02)) or 7r 2 (i(ai)) ^ 7r 2 (i(a 2 
where 7Ti and 7r 2 are projections over first and second component respectively. 
We suppose the first case (the second is symmetric) and so 7Ti o i : Sa — > Sb is 
injective. Moreover, since 

7Ti O L O G_A = 7iT O GgxC 1 — Gb O 7Tl O i, 

we conclude that A B. 

Second, we show that A < s B x C implies either A =^ s B or A =^ s C which is 
sufficient to prove the theorem for =^ s . Let I a be the set of states that can be 
reached after one step of A (formally, Ia = Ja{Sai ■ ■ ■ 1 Sa)) and Is and Iq be 
similar sets for B and C. 

• We first suppose that B and C are such that each uniform configuration 
is either a fixed-point or without any uniform antecedent. If ai I a then 
for any b and c such that s(b, c) = ai we have either b $ I® or c ^ Iq- We 
suppose the first case (the second is analogous) and then we have A <q B 
where C : Sb — >• Sa is defined by ((b) — a\ and Q(x) = a 2 for i/t. 

If a 2 ^ I a we apply the same reasoning and so we are left with the case 
I A = Sa- Since pairs of Sb x Sc arc 2-colored via s, there must be 
two pairs of different colors which agree on a component. Suppose it 
is the first component (the other case is symmetric), we have 61, & 2 £ Sb 
and c G Sc such that s(&i,c) = ai and s(& 2 ,c) = a 2 . Consider the set 
X = {(6i,c),(6 2 ,c)}. Since s o (Gg x G c ) = Ga s and s(A z ) = S\ and 
I A = Sa we necessarily have s(Sb, d) = Sa where d defined by d = Gc(c). 
d is quiescent by hypothesis on C. So we have A <q B with ( : Sb ^ Sa 
defined by £(x) = s(x, d) (( is onto by choice of d). 

• Now suppose that the hypothesis on B and C are not fulfilled. Then, if 
t = \Sb\\ x \Sc\l, both £>* and C* are guarantied to fulfill the required hy- 
pothesis (because any uniform configuration is either in a cycle of uniform 
configurations, or without uniform antecedent arbitrarily far in the past). 
Since A 1 < (B x C) = B* x C*, it suffices to apply the previous reasoning 
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on A 1 , £>' and C* to conclude either A 1 <j £>* or .4* < C*. In either case 
the theorem follows. ■ 



The Cartesian product operation is not a supremum in any of the quasi- 
order. In fact these quasi-orders don't admit any supremum or infinimum op- 
eration as shown by the theorem below. Recall that an upper semi-lattice is a 
partial order structure < equipped with a 'sup' operation such that: 

a < x and b < x sup(a, b) < x 

The definition for lower semi-lattice is dual (replace 'sup' by 'inf and any rela- 
tion x < y by y < x). 

Theorem 4.2. Let =^ be a relation among =^! s and =4 m . Then the ordered 
structure [AC/ ~, =^J is neither an upper semi-lattice, nor a lower semi-lattice. 

Proof. Let A2, A3, »42,3 and A24 be products of translations with characteris- 
tic sequences (2), (3), (2,3) and (2,4) respectively. Theorem 13.41 and lemma l3~4l 
show that they induce the following structure in =^ : 




where an arrow from A to B means B =^ A and if B =4 C =^ A then either B ~ C 
or C ~ A. This shows that the pair Ai, A3 has no supremum and that the pair 
^2,3) ^2,4 has no infimum. ■ 

4-2. Ideals and Filters 

Although the structures (AC, =4) studied in this paper are not semi- lattices 
(see above), many classical properties of cellular automata are nicely captured 
through ideals and filters. Well-known in lattice theory and algebra, the notions 
of ideal and filter can also by defined for an arbitrary (quasi-)ordered structure 
Q. For the structure (AC, =4), an ideal I is a set of CA such that: 

• if A 6 I and B 4 A then Bel; 

• for any A, B E I there is some C 6 I such that A 4 C and B 4 C. 

Moreover, / is said principal if there is some Ai such that A E I A =^ Ai- 

The notion of filter and principal filter are dual of ideal and principal ideal 
(replacing all X 4 Y by Y 4 X). 

Given a set I of CA, the three following conditions are sufficient for / to be 
an ideal for the simulation =4i (resp. =^ s , or =^! m ): 

1. A E / •<=>■ A^ E I for any transform a, 

2. if B E I and A C B (resp. A < B, or A <C B) then A el, 

3. if A E I and B E I then Ax Be I. 

Most of the proofs below follow this scheme. 
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4-.2.1. Dynamical properties 

The following theorem shows that several dynamical properties of global 
rules of CA correspond to ideals in the quasi-orders. A CA is nilpotent over 
periodic configurations if there exists a spatially periodic configuration Co such 
that all spatially periodic configurations lead in finite time to cq. 

Theorem 4.3. Let ^ be a simulation relation among =4i, and =4 m . The 
following sets of CA form ideals of (AC, ■' 

• surjective CA, 

• reversible CA, 

• CA which are nilpotent over periodic configurations. 

Proof. First, from the point of view of global maps, a geometric transform con- 
sists in iterating or composing with bijective maps. So the properties of being 
surjective or reversible are left unchanged by geometrical transforms. Besides, 
geometrical transforms map periodic configurations to periodic configurations, 
cycles of configurations to cycles of configurations (possibly reduced to a sin- 
gle configuration), and attraction basins of such cycles to attraction basins of 
cycles. Hence, nilpotency over periodic configurations, which is equivalent to 
the existence of a temporal cycle having all periodic configurations in its at- 
traction basin, is preserved by geometrical transforms. By similar reasoning on 
the phase space, it is straightforward to check that A is nilpotent over periodic 
configurations if B is and A Q B or A < B. And Ax Bis nilpotent over periodic 
configurations if both A and B are. So nilpotency over periodic configurations 
induces an ideal for 

It is also clear that surjectivity and reversibility are preserved by carte- 
sian product. Now suppose A < n B. If B is surjective then so is A since 
Ga o 7f = 7f o Gg and 7r is by definition surjective. If B is reversible, consider any 
map (f) such that 7r o cf> = Id and let A-i be the CA over state set Sa and defined 
by the global map G — W o Gg 1 o <j> (it is a shift-commuting continuous map). 
Since Ga ott = tt o Gb, one can check that Ga ° G = Id so A is reversible. 

Finally, suppose A C t B. If B is reversible then A is also reversible since 
I o Ga = Gb ° T and i is by definition injective. If B is surjective, then so is A 
because B being injective over finite configurations (Moore-Myhill theorerrjfl) A 
is also injective over finite configurations (i maps finite configurations to finite 
configurations ) . ■ 

Theorem 4.4. Let A and B be two reversible CA and =4 be a simulation relation 
among =^ s and =^ m . If A =^ B then A^ 1 =4 B . 

Proof. First, it is straightforward to check that the inverse of geometrically 
transformed instances of A are transformed instances of the inverse of A. Using 
what was shown above concerning reversibility, it is thus sufficient to prove the 
two following properties: 



3 In [T3l . one can find the following theorem: a CA is surjective if and only if there is no 
pair of finite configurations (i.e. uniform except on a finite reg ion) having the same image. 
The original formulation of the Moore-Myhill theorem [5f|, [2a| supposes the existence of a 
quiescent state. 
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• A C t B implies y^ 1 C 4 B~\ 

• i< 9 6 implies „4 _1 < g 
In the first case we have: 

Gb oT = To Ga => T = Gg 1 o I o Gj\ => I o G^ 1 = Gg 1 o l 

each equality being true on S\. In the second case we have: 

Ga°9^9°Gb^Ga°9° Gg 1 = g^go Gg 1 = G^ 1 o g 

each equality being true on Sg. ■ 

One immediate consequence of the theorem is that if two reversible CA 
are equivalent then their inverse CA are also equivalent. What is not obvious 
however is whether the inverse CA are necessarily in the same class as the initial 
CA. 

Open Problem 1. Consider any simulation relation and ~ the associated equiv- 
alence relation. What are the reversible CA F such that F ~ F^ 1 ? 

At the time of writing we have no example of a reversible F with F ^ F . 

Theorem 4.5. Let =4 be =^ or =<! TO . Then the ideal of reversible CA is principal: 
there is a reversible CA A such that 

B reversible •<=>• B =^ A. 

Proof. In |9], a reversible CA B able to simulate any reversible CA is con- 
structed. The notion of simulation used is included in and therefore in =4 m . 
The implication => is thus proven and the converse implication is proven by 
theorem 14.31 ■ 

For the ideal of surjective CA, the principality is still an open problem in 
dimension 1. 

Open Problem 2. Is the ideal of surjective CA principal, and for which sim- 
ulation quasi-order? 

Limit sets of C A have received a lot of attention in the literature 0, [IH, [Tljj . 
The limit set of A is the set fix of configurations having predecessors arbitrarily 
far in the past, formally: 

n A = f)G A (s A ). 

t 

The next theorem shows that the class of CA with a sofic limit set is nicely 
captured by =^ s . 

Theorem 4.6. The set of CA with a sofic limit set is an ideal for =<! s . 
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Proof. For CA of dimension 1, having a sofic limit set is equivalent to having a 
regular limit language [38] . It is clear that this latter property is left unchanged 
by geometrical transforms (the limit language is not affected by iterations and 
shifts, the regularity of the language is not affected by packing). Hence, it 
is sufficient to show that if B has a regular limit language and A < g B then 
A also has a regular limit language. Since regular lang uag es are closed under 
substitution (a classical result which can be found in [14]). it is sufficient to 
prove that = This last assertion is a direct consequence of A < g B, 

since the following equality holds by recurrence on t: 

g{G t B {Sl))=G t A {S\). 



Open Problem 3. Let =<! be =^ or =4 m . Is there a ^-universal CA with a sofic 
limit set? 

4-2.2. Topological dynamics 

The properties considered above are purely dynamic: they can be expressed 
as structural properties of the phase space with the reachability relation only. 
We now consider properties from topological dynamics: they are expressed with 
both the reachability relation and the topology (Cantor distance) of the space 
of configurations. We will show that many of them correspond to ideals of the 
simulation quasi-orders. 

The properties we will consider are derived from the equicontinuity classifi- 
cation of P. Kurka [19]. Let A be any CA with state set Q and global rule G 
and denote by d the Cantor distance over Q z . 

• x € Q z is an equicontinuity point for A if 

Ve,3(5,Vy g Q z : d(x,y) < S ^\/t,d(G t {x),G t (y)) < e. 

• A is sensitive to initial conditions if 

3e,\/6,Vx e Q z 3y 6 Q z 3t : d(x,y) < 5 and d(G t (x),G t (y)) > e. 

• ^4. is (positively) expansive if 

3e,Wx,yeQ z :x = y yt,d(G t (x),G t (y)) < e. 

The classification of P. Kurka is the following: 
K\ is the set of CA for which all configurations are equicontinuity points, 
K2 is the set of CA having equicontinuity points, 
K3 is the set of CA sensitive to initial conditions, 
K4 is the set of expansive CA. 
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The weakness of this classification is its lack of shift-invariance: the identity 
and the elementary translation belong to different classes (Ki and K3 respec- 
tively) . Several attempts have been made to overcome this problem by changing 
the topology @. More recently, a new approach has been proposed [j^: the 
Cantor topology is conserved (with all its good properties) but the topological 
properties are enriched with a new parameter (a velocity) which is used as the 
reference direction of information propagation in space-time. The original defi- 
nitions of P. Kurka are thus obtained by choosing velocity 0, but now identity 
and elementary translations are assigned to the same class (with different veloci- 
ties). This directional dynamic approach is more suitable for our study since, by 
definition, the equivalence classes of any of our quasi-orders are shift-invariant. 
We will define 4 classes based on the existence of some direction for which some 
dynamical behavior is observed. 

We say that A is a reseating of B if there are transforms a and /3 such that 
j[i a ) = _ w e then consider the following 4 classes: 

• the set 71 of CA which are a rescaling of some equicontinuous CA, 

• the set 72 of CA which are a rescaling of some CA having equicontinuity 
points, 

• the set 79 of C A which are not in 72 , i- e. C A which are sensitive in every 
direction^], 

• the set T4 of CA which are a rescaling of some (positively) expansive CA. 




Figure 2: Four kinds of topological dynamics. 

Figure [2] is justified by the following theorem. 
Theorem 4.7. We have the following inclusions: 

1. 7i c 75, 

2. % c %. 

Moreover, each of the sets T\, 7i \ 7i, 7jj \ 7i and 71 is non-empty. 

Proof. The first inclusion follows from definitions. The second follows from 
proposition 3.2 of [32| which asserts that the set of directions with equicontinuity 
points and the set of expansive directions cannot be simultaneously empty. 



4 For one-dimensional CA, the set of sensitive CA is the complement of the set of CA having 
equicontinuity points (see HHP . In [32t , this complementarity is shown for any direction. 



20 



Non-emptyness of 71 and 7a follows from the existence of equicontinuous 
(e.g. the identity) and (positively) expansive CA (e.g. Z^). Moreover, any CA 
having an equicontinuity point which is not equicontinuous (e.g. the CA of local 
rule <5 max (a, b, c) — max(a, 6, c)) cannot be in 7i (equicontinuity is preserved by 
rescaling), so it is in 7i \ T\. Finally, a\ x cr_i 6 75 \ %■ Indeed, for <r x cr _1 , 
any direction is either a direction of right-expansivity or a direction of left- 
expansivity, neither both. So a x ct _1 ^ Ta- Finally, a x er -1 ^ Ti since, 
by proposition 3.2 of [32j . no direction of (left or right) expansivity can be a 
direction with equicontinuity points. I 

Theorem 4.8. 

1. 71 is an ideal for any simulation =^ among =<, S7 =^ and =^ m ; 

2. T2 is an ideal for =^ s ; 

3. 71 is an ideal for =<^. 

PROOF. First, consider any A, B £ T\. Then there are CA A' and B' which 
are both equicontinuous and ^-equivalent to A and B respectively. Then, if 
C = A' x B' we have C 6 T\ and by theorem [273] we have both A^C and B ^C. 
The same reasoning can be applied to 71 and T2 ■ Thus we have shown the second 
condition of the definition of ideals for the three properties considered here. 

To conclude the theorem, and since the three properties considered are by 
definition invariant by rescaling, it is sufficient to prove: 

• if A C B or A <1 B then B equicontinuous => A equicontinuous; 

• if A <! B then B has equicontinuous points A has equicontinuous points; 

• if A C B then 23 expansive A expansive. 

The first assertion follows from the characterisation of equicontinuous CA 
as ultimately periodic CA [Tjjj]. 

For the second assertion, if A < T B then for all x, y £ S*| we have the inequal- 
ity d(ir(x), ir(y)) < d(x 7 y). Moreover, for any yi £ S\ there is some yi £ Sg 
such that 7r(?/2) = yi and d(ir(x), j/i) = d(x, 1/2) (choose 2/2 so that it equals x on 
the cells around position 0). Hence, if x is an equicontinuous point for B then 
ir(x) is an equicontinuous point for A. 

Finally, for the third assertion, it is sufficient to notice that the property 
of expansivity is defined by a formula using only universal quantifications on 
configurations so it remains true on a subset of configurations. I 

72 is not an ideal for =<;,; and neither for =4 m as shown by the following 
example. 

Example 4.1. Consider B £ T3 of radius 1 and let A be the CA with radius 1, 
states set Sa — SbU {M} (with M ^ Sb) with local rule f^ defined by 



B C A so B =4i A. However A £ T2 since the configuration u M w is an equicon- 




tinuous point. 
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Notice also that 7ei cannot be an ideal because a\ x cr_i £75 simulates 
o- E Ti. 

Open Problem 4. Are there A ^ T4 and B £ 7i smc/i </ia< .4 =^,5 £> (i.e. £/ie 
simulator CA is expansive up to reseating but the simulated CA is not expansive, 
even up to reseating)? 

4-3. (Un) decidability 

The fact that many properties related to the simulation quasi-orders are 
undecidable is no surprise. For instance the nilpotency property, which is an 
undecidable problem [18| , corresponds to an equivalence class in the three quasi- 
orders (theorem [XT]). However, there are non-trivial properties of these quasi- 
orders which are decidable (see below) and the edge between decidable and 
undecidable properties is hard to catch. 

In this section, we consider two kinds of problems in simulation quasi-orders: 
lower bounds (being above some fixed CA or set of CA) and upper bounds (being 
simulated by some fixed CA or some CA from a fixed set). 



Theorem 4.9 ([23]). The set of CA of radius 1 with a spreading state and 
nilpotent over periodic configurations is not co-recursively enumerable. 

Theorem 4.10. Let A be any CA which is not nilpotent over periodic config- 
urations. Let =3! be either =<!j or =4 m . Then the set of CA B such that A =4 B is 
not co-recursively enumerable. 

Proof. We describe a computable transformation which, given a CA C of radius 
1 with a spreading state, produces a CA B with the following properties: 

• if C is not nilpotent over periodic configurations then A =4 B; 

• if C is nilpotent over periodic configurations then so is B. 

The theorem follows by theorem l4.3l since we have reduced the problem 'A =4 BT 
to the problem of nilpotency over periodic configurations (reduced to CA of 
radius 1 with a spreading state). 

We now describe the construction of B from C. Suppose C has a spreading 
state q. B is the CA of radius 1 and states set Sjg = (Sc \ {q}) xS^U {q} with 
local rule fs defined by: 

{It! u \t< h \\ . r ja,b,ce S B \{q} and 
[fe{ai,bi,ci) ?= q, 
q in any other case, 

where a,, bi and Ci represent component i of a, 6 and c. Any periodic configura- 
tion c of B either leads to the uniform configuration g, or contains a periodic con- 
figuration of C in its first component. Hence, if C is nilpotent over periodic con- 
figurations, then so is B (because q is precisely the spreading state of C). If C is 
not nilpotent over periodic configurations, then there is a word u 6 (Sc \ {q}) m 
and an integer t > 1 such that the periodic configuration c of period u verifies 
G^c) = c. Therefore, by definition of B, we have A {mXt - Q} B^.bt.o) w hcre 
i : -> Sb is defined by: 

i(ai,a 2 , . . . , a m ) = ((01,1*1), . . . , (a m ,u m )). 
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This result shows that it is generally undecidable to know whether a CA is 
lower-bounded by a given (fixed) one. However, there are noticeable exceptions 
in one dimension for =^ s and =^ m . 

Theorem 4.11. Let =<; be either =<! s or =<! m and let A be a nilpotent CA. Then 
the problem of determining if a given B is above A for =^ is decidable. 

Proof. We are going to show that A =4 B if and only if B is not surjective 
and the theorem follows by decidability of surjectivity in one dimension [lj. 
First, by theorem 14.31 if B is surjective then A^B. Suppose now that B is not 
surjective, i.e. that B possesses some Eden word u <E for some length m. 
Then, denoting by C the CA over states set Sc = {0, 1} which is constant equal 
to 0, we have C <^ B 1 -™' 1 ^ if tt : S™ -> S A verifies tt(iu) = if and only if 
w = u. We deduce by theorem 13. II that A =4 B. ■ 

Open Problem 5. Is there a non- surjective CA A which cannot injectively 
simulate any nilpotent CA ? Is the problem of being above the class of nilpotent 
CA for injective simulation a decidable problem? 

Concerning upper-bound problems, the edge between decidability and un- 
decidability is also non-trivial. For instance, theorem 14 . 51 shows the existence of 
a CA A such that the upper-bound decision problem '£> =<! AT is decidable in 
dimension 1. 

5. Tops of the Orders 

In this section we study maximal elements of the quasi-orders. These CA 
are able to simulate any other CA 

Definition 5.1. Let =<! be any relation among =^,; 7 =^! s and =4 m . A CA A is said 
^-universal if for any B we have B =4 A. It is strongly ^-universal if it strongly 
^.-simulates any other CA. 

The notion of strong ^-universality above is exactly the same notion as 
intrinsic universality defined in section 5 of Q and has already been considered 
several times in the literature (see (3ll | for a survey). In fact, strong and general 
universality are the same notion for =<!j and =4 m . 

Theorem 5.1. There exist strongly ^i-universal CA and all =4i-universal CA 
are strongly ^i-universal. The same is true for =^ m . 

Proof. For the existence of strongly ^-universal CA, see [3l|. The theorem 
follows by application of theorem 12 of [8|. ■ 

Of course, any ^-universal is also =<; m -universal. The converse is an open 
problem. 

Open Problem 6. Do the notions of ^i-universality and ^ m -universality co- 
incide? 

Concerning =<t s , the situation is different: no CA is strongly =^ s -universal 



5 The proof of this fact was suggested by G. Richard. 
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Theorem 5.2. There is no strongly =^ s -universal CA. 

Proof. Suppose for the sake of contradiction that there is some strongly =^ s - 
universal A. Consider a uniform configuration c of A. There is n such that the 
orbit of c under A contains n different configurations (the orbit is ultimately 
periodic) . Now consider B with n + 1 states such that its uniform configurations 
are all in the same cycle of length n + 1. By hypothesis, for any B there is some 
geometric transform a such that B < s A^ . Let d be the corresponding config- 
uration of c for A^ a ' . The orbit of d contains at most n different configurations 
and it is therefore the same for the orbit of s(d) under B. But s(d) is necessarily 
uniform and we get a contradiction with the choice of B. I 

The theorem 12 of don't apply for =^! s . However, we are not able either 
to construct a ^-universal CA, nor to prove that there is none. 

Open Problem 7. Is there a =4 s -universal CA? 

For the rest of this section, we consider only =4i and =^ m . 

5.1. On the Way to the Top 

Universal CA are not hard to construct and the property of being universal is 
recursively enumerable since simulation relations considered here are recursively 
enumerable. However universality is not co-recursively enumerable as shown by 
the following theorem. The case of ^-universality was proven in [30| . Using 
theorem 14.101 the proof below is direct and includes the case of =^ m . 

Theorem 5.3. The set of ^.{-universal CA is not co-recursively enumerable 
and neither is the set of =4 m -universal CA. 

Proof. There exists a CA which is ^-universal but not nilpotent over periodic 
configurations. To see this consider any universal CA and add a new state 
which is spreading: the resulting CA, say A, contains at least two disjoint 
periodic orbits of periodic configurations and is thus not nilpotent over periodic 
configurations. The theorem follows by application of theorem 14.101 to A since 
A is by construction both ^-universal and =<; m -universal. ■ 

This result has some consequences on the structure of simulation quasi-orders 
'near' the top. The following theorem shows that a non-universal CA is always 
'infinitely far' from the class of universal ones. 

Theorem 5.4. Let =^ be =^.; or =^ m . And let ii be the set of ^-universal CA. 
Then we have: 

1. AxBeii <s=> AeU orBeii, 

2. ifA^U. then there is B <j£ ii with A ^ B but B^A. 

Proof. 

1 . By theorem 12.31 we have A =^ A x B and B =4 A x B which proves one di- 
rection. Moreover, there exists C 6 ii with 2 states only 0, [28|. If we 
suppose Ax B € il then, by theorem 15.11 it strongly simulates C. Hence, 
by theorem 14.11 we have either C =4 A or C =4 & and thus either A G it or 
Be ii. 

2. Let A ^ ii. If A was such that C =^ A for all C il then the complement of 
il would be the set {C : C =^ .4} and il would be co-recursively enumerable 
contradicting theorem 15.31 So there is C il with C^A. To conclude the 
proof it is sufficient to choose B = Ax C (theorem I2.3[) . ■ 
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5.2. Necessary But Not Sufficient Conditions 

The purpose of this section is twofold. It compares the notions of universality 
defined above to other definitions of the literature and, by doing this, presents 
tools and techniques to prove non-universality of some CA (other proofs of non- 
universality for other purposes are developped in section [5]). 

One of the techniques we use to ensure that some CA is not universal yet 
achieving some behavior B, is to add a spreading state and let the CA generates 
this state if it detects somewhere that the current configuration doesn't corre- 
spond to a 'legal' configuration, i.e. a configuration occuring normally when 
producing the behavior B. Proofs of non- universality with this technique rely 
on the lemma below. Before stating and proving the lemma, we need to give 
some precisions on spreading states and sets of configurations 'supporting' a 
simulation. 

First, the notion of spreading state is sensitive to the choice of the syntactical 
representation of the CA because it depends on the choice of the neihborhood. 
In the sequel we say a CA A has a spreading state k if any cell changes to state 
k when k appears in its minimal neighborhood (i.e. the minimal set of cells 
upon which the local rule effectively depends). 

Second, given a relation of the form A C< B^ m ' l,t ' z \ there is an isomorphism 
between {A,S%) and (B {m ' 1 ^' z \(i(S A )) 1 ') as dynamical systems. At the level 
of B, the configurations involved in this relation is the set X of configurations 
made of infinite concatenation of elements of i(S A ) C Sg l (viewed as words of 
length m over alphabet Sg). This kind of sets are called block- sub shifts and 
discussed in more details in section 3.2 of Q. In the sequel, such a set X is 
called the support of the simulation. 

Lemma 5.1. Let A be a CA without spreading state and B be a CA with a 
spreading state K. If B strongly ^4 rn -simulates A, then the support X of the 
simulation cannot contain k. 

Proof. By hypothesis, there are parameters m, t, r, z and a CA C such that 

A<„CQi B< m ' T >*-*>. 

By choice of B, B ( ~ m ' T ^ admits n m as spreading state. Moreover, by definition of 
<C, the minimal neighborhood of A is included in the minimal neighborhood of 
g(m,T,t) _ Thus, if k appears in some configuration of X then the state 7r(i _1 (K m )) 
is a spreading state for A because K n also appears in X for arbitrarily large 
n. ■ 

We first study how embeddings of Turing machines into CA can relate the 
notions of universality for Turing machines to the notions of universality derived 
from quasi-orders as defined above. 

An embedding of a Turing machine A4 into a CA A is an embedding of the 
instantaneous descriptions of Ai into configurations of A such that instanta- 
neous descriptions of successive steps of M. corresponds to successive steps of 
A via the embedding. We don't give any formal definition of embedding since 
we will never prove negative results (i.e. assertions of the form 'there is no 
embedding of M. such that...'). However, the embeddin gs w e use in the sequel 
are classical and already appeared in the literature (see [34J). 
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Theorem 5.5. For any Turing machine M, there exists a CA A which embeds 
M but is not ^4 m -universal. 

Proof. Let M = {Sm,Qm,4>m) where Sm is the set of states of M, Qm is 
the tape alphabet, and 

4>m '■ Sm x Qm -> Sm x Qm x {— 1, 0, 1} 

is the transition function of M.. We construct a CA A over state set 

Sa = Qm x {<-, U Qm xSjmU {k} 

where —> and «— are states not already in Sm ■ Each cell of A corresponds to 
a tape position of M: it contains a letter from the tape alphabet and either 
a head with its current state or no head but an indication or — > telling in 
which direction to find the head. On configurations containing a single head, 
A mimics transitions of M. step by step as expected. Thus, A embeds M.. In 
addition, A checks that <— a never occur to the left of a state from Sm ora4 
(and symmetrically for — >). If the check fails, then the state n is generated and 
spreads. 

This construction ensures that, for any initial configuration c, if the orbit of 
c never contains an occurrence of k then it contains at most one head. Hence, 
these orbits are such that at any time step state changes occur on the neigh- 
borhood of at most one position (a head move involves a state change in two 
adjacent cells). 

Now suppose that A is =<; m -universal and consider the CA B = a\ x er_i. A 
strongly simulates B (theorem 15. ip . Since B has no spreading state, then the 
set X of configurations of A on which the simulation occurs never contains k. 
We deduce that all orbits of configurations from X have the property described 
above. This implied that B is such that on all its orbits, at most two cells change 
their states between two steps: this in contradiction with the choice of B. ■ 

Turing-universality of cellular automata is a fairly vague notion in the liter- 
ature. We don't give a formal definition here since we won't prove any negative 
result concerning Turing-universality. We just consider that a CA able to embed 
a universal Turing machin<|f| is Turing-universal. 

We can choose M to be universal in the previous theorem (theorem [53]). In 
this case, since the embedding used in the proof ensures that M. is simulated in 
real time by A, we deduce that the following problem is P-complete: 

Input: a state q e Sa, an integer t > 1, and a word u 6 S^[ t+1 where r is the 
radius of A; 

Query: do we have ^4*(u) = ql 

This problem of finite triangle computation has been considered several times in 
the literature and it has been proven that it was P-complete for particular CA 



12|>|27( . This notion of complexity inherited from sequential computation theory 



fails to capture the notion of universality associated to simulation quasi-orders. 

Corollary 5.1. There exists a CA which is Turing-universal and P-complete 
but not =4 m -universal. 



6 We don't give any formal notion of universality for Turing machine either. In fact, we 
only need to suppose the existence of at least one universal Turing machine. 
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6. Induced Orders 

This section aims at studying particular CA or sets of CA for the ordered 
structure they induce in the simulation quasi-orders. While studying various 
properties of the quasi-orders in the previous sections, we have already estab- 
lished the existence of several induced infinite structures. 

For instance, theorem 15.41 allows to construct an infinite strictly increasing 
chain of non-universal CA starting from any non-universal CA for the quasi- 
orders associated to =4i and =<; m . Besides, theorem 13.41 implies the existence of 
inifinite chains in the three quasi-orders studied in this paper. 

Section 16.11 below gives a way to construct chains of length u) + lo and an 
hint about the existence of chains of length u> x lo. However, we leave open the 
question of the longest chain induced in any of the quasi-orders. We don't even 
know if one of them admits a dense chain. 

Open Problem 8. Does one of the quasi-orders admit a dense induced order? 

6.1. Limit Cartesian Product 

We have seen in theorem 15.41 that if A is not universal, then A x A cannot 
be universal. Therefore, no finite Cartesian product of A with itself can be 
universal. Therefore, we have a chain of non-universal CA: 

.4=^.4x.4^.4x.4x.4^--- 

For some A, the chain collapses in a single equivalence class, e.g. if A is a 
translation (see lemma I3.4j) . However, the following theorem shows that for 
some A, the chain is strictly increasing. Moreover, A can be chosen so that it 
embeds any Turing machine. 

Theorem 6.1. For any Turing machine M., there is a CA A which embeds M. 
and such that for any 1 < n < m, one has: 

A x ■ ■ ■ x A Ax ■ ■ ■ x A . 

m n 

Proof. Le A be the CA constructed in the proof of theorem 15.51 We can 
suppose that M. is such that it can produce infinite sequences of left move of 
its head when started from a special state (not the initial state) over a blank 
tape, and more precisely that the sequence of moves leaves the tape blank. If 
Ai does not have this property, just add some states to achieve this behavior. 
We can suppose the same for right moves. 

Denote by B rn the product of m copies of A and by B n the product of 
n copies. Suppose for the sake of contradiction that B m =4 m B n q. We can 
construct for any set of positions z\, . . . , z m a configuration c of B m such that 
for all i the ith component contains a correct instantaneous description of M 
where the head is at position Zi in a state suitable to generate an infinite sequence 
of left or right moves (as supposed above). Now let d be a configuration of B n 
corresponding to c via simulation. First, if some component i of c' contains a 
spreading state, it will spread and, after some time t, will be present at some 
position where the configuration (c) contains no head, but only a blank 
tape symbol on each component. This means that blocks of blank tape symbols 
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in B m can be simulated by blocks of B n where the ith component is a block 
of spreading states. Considering again the orbit of c, we deduce that it can 
be simulated by a configuration c" where the ith component is everywhere a 
spreading state except at a finite number of positions. Thus after some time, 
the ith component will become uniform and constant. It is then straightforward 
to show that it is useless for the simulation and that in fact B m on c can be 
simulated by only n — 1 copies of A. 

Applying the reasoning inductively, we can therefore suppose that no spread- 
ing state appears on any component in the orbit of the configuration c' defined 
above. Since, the orbit of c is such that there are m distant positions where some 
states change at each step, it must be the case in the orbit of c'. Since, n < m, 
there must be some component with two heads and therefore a spreading state 
must appear after the first step: this is in contradiction with what we have just 
supposed. ■ 

For the CA A of the previous theorem, we can ask if the infinite chain of 
Cartesian products is upper-bounded by some non-universal CA, or if any CA 
able to simulate each product of the chain is necessarily universal. One can 
imagine that for a sufficiently simple A, there is some room above the chain of 
products of A and below the class of universal CA. 

The rest of this section is devoted to the proof of a stronger result: for any 
A, there is a CA B which is able to simulate any finite product of A and such 
that B is universal if and only if A is universal. Moreover, B can be obtained 
from A constructively. Because it extends property of Cartesian product given 
by theorem 15.41 this construction will be called limit product in the sequel. If 
A is a CA, its limit product is denoted by .Aoo- 

Note: In the rest of this section we only consider the simulation =^ m . 

Without loss of generality, we can suppose that A has radius 1 ftheorem l2.3l) . 
To be able to simulate the product B of n copies of A, Aoo is made of three layers 
(its state set is a Cartesian product union a single state, which is a spreading 
state as explain hereafter): 

1. the state layer, 

2. the transport layer, and 

3. the synchronisation layer. 

It proceeds as follows. 

State layer. Each component of a cell of B is simulated by a block of three 
adjacent cells in the state layer of Aoo- More precisely, component i 
(0 < i < n — 1) of cell z of B is simulated by the block of three cells of 
Aoo beginning at position 3(nz + i). This block is referred to as B Zj i in 
the sequel. In B^, the center cell stores the ith component of the cell z 
of B and the two other are used to store temporarily the ith components 
of cell 2—1 and z + 1 . 

Transport layer. The role of the transport layer is precisely to bring states of 
ith components corresponding to cell z — 1 and z + 1 of A to the dedicated 
cells of Aoo in &z,i- Then, the transition f^(x z -i, x z , x z +i) of the ith 
component of B can be simulated locally by Aoo in B 2 ,j. Transport is done 
in parallel for any i and any z. To do this, the transport layer is made 
of a succession of particles (one every three cells), each one being able to 
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carry a state of A. Initially aligned with the center of blocks, the particles 
move in parallel according to a cycle of five steps: 

1. move right by 3n cells and read the state seen on the state layer; 

2. move left by 3n — 1 cells and write the memorized state on the state 
layer; 

3. move left by 3n + 1 cells and read the state seen on the state layer; 

4. move left by 3n — 1 cells and write the memorized state on the state 
layer; 

5. move 1 cell right and apply local rule f a on state layer at the current 
position; 

Synchronisation layer. The role of the synchronisation layer is to orchestrate 
the cycle of particle moves and it must be able to do it for arbitrary 
large values of n (simulating arbitrarily large Cartesian products of A is 
sufficient to simulate all products of .4). It contains a flag that can take one 
of the four indications 'left', 'right', 'read', 'write' and 'transition'. The 
flag is changed everywhere synchronously according to a cycle suitable to 
ensure that particles of the transport layer produce the cycle described 
above when they follow the instruction given by the flag. 

We now describe in detail the synchronisation layer. Denote by u n the flag 
sequence mentionned above in the simulation of a product of n copies of A, and 
let E be the set of flag states. 

Theorem 6.2. There is a CA C with a spreading state k and a map tt : Sc — > E 
such that C is not ^ m -universal, and, for any configuration c G S®, one of the 
following property is true: 

Cycle: at each time in the orbit of c, all cells have the same image by n and 
the sequence with time of this common image is periodic of period u n for 
some n; 

Frozen: at each time in the orbit of c, all cells have the same image by n, but 
this common image remains constant after a certain time; 

Error: the spreading state appears at some time in the orbit of c. 

Moreover C is such that there are configurations having the 'cycle ' property above 
producing period u n for arbitrarily large n. 

Proof. First, notice that flag changes in the sequence u n are separated by a 
number of steps which is either constant (independant of n), or of the form 
3n + c with c a constant (we can suppose c > without loss of generality) . To 
simplify notations, we will suppose in this proof that u n alternates between two 
values and 1 every 3n steps. Adapting the proof for the real u n is just a matter 
of adding a finite set of special states to deal with constants. 

The proof is based on a reversible solution B to the firing squad synchro- 
nization problem proposed by K. Imai and K. Morita: in [l6j |. they construct a 
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reversible CA B with a subset of states F (the firing states) such that for any 
n, there is a periodic configuration c„ verifying: 




\c n ) e F z 

(c„) e (S B \ Ff for all t, < t < 3ra. 



• Gt 



B 



Without loss of generality, we can suppose that B and its inverse are syn- 
tactically represented with the same radius. We now define a CA Co of radius 
r, with states set Se = Sjg x Sg x {0, 1}, and with transition function: 



where x( a ; b) equals 1 — b if a £ F and b else. Intuitively, on configurations 
where the third component is uniform equal to b, Co mimics B on the first com- 
ponent and £> _1 on the second one if b = 1 or the converse if b = 0. Moreover, 
the value of b is switched each time the component playing B encounters a firing 
state. Hence, if we choose for tt the projection on third component, Co started 
from configurations c„ has the property 'cycle' and produces the periodic se- 
quence u n . 

We now enrich Co with a spreading state which is produced each time one 
of the following local checking fails: 

• the third component {0, 1} must be uniform; 

• for the two first components, a state from F (firing state) must always be 
surrounded by states from F only; 

• states from F are forbidden on the second component if b — I and states 
from F in the first component are forbidden if b = 0. 

The third condition ensures that in the case of a 'cycle' regime (firing states ap- 
pearing infinitely often) , the period is equally divided between steps where b = 
and steps where 6=1. To ensure that such a 'cycle' regime always produces an 
alternance of exactly 3n zeros and 3n ones, we add a component implementing 
a counter modulo 3: the value of this component is incremented modulo 3 at 
each step (whatever the context) and a spreading state is generated if a cell 
contains a firing state and the counter is not modulo 3. Denote by C the CA 
obtained and consider any configuration c. If no spreading state appears in the 
orbit of c, then the third component is uniform. If it changes of state only a 
finite number of times, then we are in the 'frozen' regime. If there are infinitely 
many changes, it follows from the discussion above that the conditions of the 
'cyclic' regime are fulfilled. 



7 In lial . the main concern is synchronisation of finite segments of cells surrounded by a 
quiescent state. To extend the property to infinite configurations, it is crucial that "garbage" 
(which must be conserved to ensure reversibility) do no spread outside the initial segment. The 
solution of K. Imai and K. Morita has precisely this property as it is cxplicitely mentionncd 
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in [3. 
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To conclude the theorem, it remains to prove that C is not ==! TO -universal. 
Suppose for the sake of contradiction that it is and let ti be any universal CA 
without spreading state and consider the set X of configurations of C which is the 
support of the strong simulation of U (C strongly simulates U by theorem 15. ip . 
X cannot contain any occurrence of the spreading state (by lemma I5.ip . it 
implies that all configurations of X have a uniform third component. But, on 
such configuration, the dynamics of C is reversible. Hence U is reversible: this 
is a contradiction with its universality by theorem 14.31 I 

The synchronization layer of limit products is exactly the automaton C of 
the previous theorem, except that the spreading state of C now becomes a global 
spreading state. Before establishing the main result of this section, we give more 
details concerning the state layer and the transport layer of Aoo ■ 

The state layer is made from state set Sa x {L, C, R} where L, C and R 
are states to identify explicitly the role of each cell in each block B^j: C for the 
cell storing the ith component of cell z of A and L and R to temporarily store 
states of ith component of cells z — 1 and z + 1 respectively. 

The transport layer is made from state set Sa U {_!_} where _!_ is the state 
used to separate particle carrying a state from Sa- 

So the states set of Aoo is: 




state transport synchronization 



In addition to the behavior described above, Aoo does the following local 
checkings and generates the spreading state k if one of them fails: 

• the second component of transport layer must be periodic of period LCR; 

• the transport layer must contain an alternance of one state from Sa an d 
two states _L; 

• when doing read and write operations, the particles of the transport layer 
must be aligned with the right type of state in the state layer: 

— type C when reading, 

— type R when writing at step 2, 

— type L when writing at step 4; 

• when the synchronisation layer says 'transition', check that the particules 
are aligned with cell of type C in the state layer. 

All those checkings ensure the following property: if no spreading state is 
generated and if the component layer produces a correct cycle of instructions, 
then the behavior of the state layer is equivalent to the behavior of some Carte- 
sian product of A (up-to some rescaling) . 

Before stating the main theorem, we establish a simple yet useful lemma 
saying that if A simulates B with support X, then everything A can simulate 
using a support included in X can also be simulated by B. 

Lemma 6.1. Let =^ be either or =4 m . Let A and B be such that the simu- 
lation A =4 B occurs on a support X of configurations of B. If B ^-simulates C 
on a support included in X , then A ^.-simulates C. 
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PROOF. We consider the case where is By hypothesis, we have A^ a ' Q B^ 1 ^ 
on support X and Qj B^ 2 ' on support FCI. Now, let m a , m^, mp 2 and 
rrij be the packing parameters of transforms a, /3i, fa and 7 respectively. The 
injective maps i and j induce two injective maps ip 2 and jp t with the following 
domains and ranges: 



l/3 2 . D A —r^B 



Therefore = i^ 1 o is a well-defined injective map from s™~* m Pi into S^°' m?2 _ 
Now define the transforms ?y a and ?y c to be the composition of a and $2 , and of 
7 and /3i respectively. Then we have C^ 1 "' ^l' 7 '"). 

The extension of the previous reasoning to =^ m is straightforward. ■ 

This lemma together with lemma 15.11 is the key to a kind of 'self-checking' 
simulation used in the construction of the limit product (and re-used in sec- 
tion A 'self-checking' simulation of B by A is standard simulation of B by 
A on some support X with the additional property that A 'checks' locally on 
any configuration that it belongs to X, and triggers some pathological behavior 
(typically a spreading state) in case of check failure. Hence any possible strong 
simulation of some C by A is such that: 

• either it has a support included in X in which case B can also simulate C 
by lemma [6TTI 

• or it must contain some c ^ X in its support in which case a spreading 
state is generated and lemma I5TT1 gives some limitation on C. 

To show that a spreading state is generated in the second case above, a 
crucial property is that the support of any simulation is by definition always 
irreducible: if u\ and 112 are words occuring in two configurations of the support, 
there exists a third configuration of the support where u\ and ui both appear 
(see section 3.2 of [§] for a more detailed discussion on supports of simulations). 

We now state the main theorem of this section. 

Theorem 6.3. For any A, its limit product Aoo is such that: 

• Ax ■ ■ ■ x A =4 m Aoo for all n > 1, 

n 

• Aoo is =4 m - universal if and only if A is =4 m -universal. 

Proof. The first assertion follows from the construction of Aoo and the detailed 
discussion above. Now suppose that „4oo is =^ m -universal and let U be any 
universal CA with no spreading state. By theorem 15.11 Aoo strongly simulates 
U:U^ Aoo {a) for some geometrical transform a. Let X denote the support of 
the simulation. By choice of 1A, the spreading state k cannot appear in any orbit 
of any configuration of X (by lemma liO]) . We deduce from theorem 16.21 that 
the synchronization component is in the same regime (either 'cycle' for a fixed 
value n or 'frozen') for all the configurations of X because otherwise, we could 
construct a configuration in X producing a spreading state (by irreducibility of 
X). 
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In the case where all configurations are in the frozen regime, the flag of the 
synchronization layer becomes constant after some time to, so the transport 
layer has the behavior of a translation (or identity) and the state layer remains 
constant, to is identical for all configurations of X (because otherwise, we 
could once more combine two configurations to produce a spreading state, by 
irreducibility of X). Then, consider a CA U + with state set Su X {0, . . . ,to} 
which has the following behavior: 

• the second component is decreased by one until it reaches 0; 

• on the first component, the local rule of U is applied, but only if the second 
component is 0. 

Since U is =<; m -universal, it can strongly simulate U + (by theorem l5TTj) : precisely, 
U+ <C l/' 7 ™' 1 '*' 2 '. Consider the set Y of configurations of U corresponding via 
simulation to the set of configurations of IA+ uniformly equal to to on the second 
component. Denote by Xy C X the corresponding set of configurations of Aoo- 
By choice of U+, we know that U simulates itself on the set of configurations 
G%° (Y). This implies that for some t' > to, Aoo can simulate U using as support 
the set of configurations (Xy)- By hypothesis, starting from such configu- 
rations, Aoo has a behavior of translation or identity on the state and transport 
layers. Since the synchronizing component evolves independently of the others, 
we deduce by lemma 16.11 that there is some CA B which is a product of trans- 
lations (corresponding to state and transport layers) such that B x C simulates 
U: this is a contradiction by theorem 15.41 since neither B (theorem [23]), nor C 
(theorem 16. 2p is universal. 

Hence, we are necessarily in the case where the synchronization layers pro- 
duce a valid cycle. Since no spreading state can be generated in the orbit of any 
configuration of X, the state layer always behaves like a Cartesian product of 
n copies of A. The value of n is in fact common to all configurations of X (as 
shown above), so we deduce by lemma l6~T1 that A x • ■ • x A simulates U and A 

n 

is therefore universal by theorem 15.41 ■ 

Of course, we can consider Aoo itself as a new candidate for taking its finite 
Cartesian products and applying the limit product construction. In fact, the 
process can be repeated forever with the guarantee that no CA ever produced 
in this chain will be universal, provided the initial CA is not. However, there 
is no reason why this infinite chain should be strictly increasing. In particular, 
even if 

.4=^.4x.4^.4x.4x.4^--- 

is a strictly increasing chain, it might be the case that Aoo is equivalent to 
Aoo x Aoo ■ Therefore we have only proven that one of the following properties 
is true: 

• there is a strictly increasing chain of length uj x lo in the quasi-order 
{AC, =4 m )i 

• for any non-universal CA A, there is a non-universal CA B such that 
A =4 m B and B x B is equivalent to B. 
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6.2. Sub-Families of Cellular Automata 

Theorem 12.31 shows that any equivalence class in any quasi-order contains 
some CA with radius 1. This fact is a direct consequence of a well-known trans- 
formation of CA with large radius into CA of smaller radius with more states 
(this transformation is called 'higher block presentation' in symbolic dynamics, 
see UK). 

It is also sometimes invoked in the literature that considering CA with two 
states only is not restrictive since there is a converse transformation that trans- 
forms a CA with many states into a CA with only two states but a larger radius. 
However, the situation is not similar to that of radius reduction since there are 
equivalence classes with no 2-states CA: e.g. Z p for any prime p ^ 2 as shown 
by theorem 13.21 and lemma 13.21 Note that the same is true for any fixed state 
set of cardinal n: the equivalence class of Z p contains no such CA provided p is 
prime and does not divide n. 

Hence, this transformation introduces a bias: the transformed CA may be 
inequivalent to the original one. Meanwhile, we know that CA with two states 
can be as powerful as CA in general since there are universal CA with two 
states only 0, [28| (for simulation relations =^ and =^ m ). More precisely, as 
we will see below, the transformation applied to a universal CA always yields 
a universal CA because the transformed CA simulates the original one. Since 
the original and the transformed CA are not always in the same equivalence 
class, one question that naturally arises is: what CA can be simulated by the 
transformed CA but not by the original one? Although it provides only partial 
answers, this section is devoted to that kind of questions, for CA with two states 
and for other families. 

Formally, given a familly J of CA, we say that a map <f> : AC — > "$ is a 
encoding of CA into familly 5 if 

VA,A<4<f>{A). 

We will only consider simulation relations and =^ m in the sequel, thus an 
encoding into $ implies that there are universal CA in A trivial example 
of such an encoding is given by $ = {U} where U is a universal CA and <j) is 
the function mapping any CA to U. We are interested in using this notion 
of encoding with families which are more 'representative' of the diversity of 
behaviors in the whole set of CA. To express this we introduce the following 
notion of faithfulness. 

Given a ^-encoding <j> : AC — > $ and a set E of CA, we say that 4> is faithful 
for E if: 

\/BeE:B^A B^ 4>{A). 

An encoding is faithful for E if the original CA and its image by the encoding 
simulate exactly the same CA in E. So, to give some evidence that a family # 
is 'representative' of CA in general, we can exhibit an encoding of CA into # 
which is faithful for a set E of CA as large as possible. When E is the whole set 
of CA, the faithfulness implies that there is a CA of family # in any equivalence 
class: this is the case for CA with radius 1. 

The next theorem gives four encodings which are faithful for U, the set of 
=^ m -universal CA. The families corresponding to these encodings were already 
defined in this paper except one: captive CA. 
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Captive CA were introduced in [35J and are denned by a simple restriction 
on the transition rule. A CA A, with state set Sa, radius r and local rule fA is 
captive if: 

Va_ r , . . . , a r £ Sa ■ jU(a-r, ■ ■ ■ ,a r ) e {a_ r , . . . , a r }. 

In the following theorem, encodings are different but their faithfulness rely on 
the same idea of 'self-checking' simulation explained above which uses lemma [5~T1 
and lemma HTT1 

Theorem 6.4. Let =<; be =^ or =<; m . For any family of CA below, there is a 
^.-encoding from CA into $ wich is faithful for the set U of =4 m -universal CA: 

• CA with two states, 

• CA in 

• CA in T~3, 

• captive CA. 

Proof. To describe the encoding for each family, we suppose A is a CA with 
state set Sa = \p>i, • ■ • , a n }, with radius r and location rule fA- 

2-states CA. Let m be an integer large enough and ip be an injective map 
from Sa to {0, l} m such that no word 4>(a) contains an occurrence of 11. 
Now define the injective map i : Sa — > {0, l} m + 4 by i(a) = 0110^i(a). Let 
r' = (r + l)(m + 4). 4>(A) is a CA of radius r' and state set {0, 1} defined 
as follows: 

• on the set X of configurations made of infinite concatenations of 
words from z(5^), <j){A) is isomorphic to A so that A C« 4>(-A); 

• everywhere else, 4>(A) generates a 1. 

The map <j> is thus an encoding of CA into 2-states CA. Now suppose 
that 4>{A) is universal and let U be a universal CA with two states and 
no spreading state which is strongly simulated by (f>(A) on support Y 
(theorem 15. 1| . If there is some y G Y with y X then 

• either there are two occurrences of 0110 in y which are not correctly 
spaced, 

• or there is a word OllOwOllO occurring in y with u i/j(Sa)- 

In any case, the image of y will contain an occurrence of 111 (because the 
above error must be seen by at least three consecutive cells) and l's will 
propagate like a spreading state which is impossible by lemma [5TT1 because 
otherwise cf>(A) ^ 3 ' 1 ' 1,0 ^ could simulate U on a support where it possesses 
the spreading state 111. So Y C X and lemma I5TT1 shows that A simulates 
U. Hence A is universal if and only if <f>(A) is. 

Captive CA. The encoding technique for ca ptiv e CA is very similar and al- 
ready appeared in a non- faithful form in [35j. Let u be the word a\ • • ■ a n> 
let # be a state not in and denote Q = Sa U {#}■ We define the 
injective map i : Sa — > Q n+3 by i(a) = ^uf^a. We then define <fi(A) in a 
way similar to the case above. Its radius is r' = (r + l)(n + 3), its state 
set is Q and its local rule is such that: 
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• on the set X of configurations made of infinite concatenations of 
words from i(Sj), <X-4) is isomorphic to A so that A 4>(A); 

• everywhere else, 4>(A) take as new state the maximum of its neighbors 
for some fixed ordering of Q such that # is the maximum. 

First, 4>(A) is captive and (f> defines an encoding of CA into captive CA. 
Second, notice that for any support of simulation Y of 4>(A), if there is 
some y GY with y g" X then, by irreducibility of Y, either there is y' G Y 
with y' g" X and y' contains a #, or # never appears in Y. In the second 
case, 4>{A) always applies a max as local rule and therefore possesses 
a spreading state when restricted to Y. In the first case, consider the 
configuration y' and z G Z such that positions z and z + 1 both see a # in 
their neighborhood and a local pattern not in X (such a z must exist by 
choice of y' and definition of X). Then </f>(„4)(y') contains the pattern 
which is spreading by definition of <fi(A). In any case we can apply the 
usual reasoning with lemma [5~T1 and lemma l6~Tl any CA without spreading 
state strongly simulated by 4>{A) is also simulated by A. So the encoding 
</> is faithful for universal CA. 

Ti and 73. For 7i, the encoding is simple: <fi(A) is just „4 with an additional 
state k which is spreading. The resulting CA <fi(A) is always in T2 since 
K 2r is a blocking word (sec 19]). Lemma [5.11 is then sufficient to prove 
that it is an encoding from AC to T3 which is faithful for universal CA. 

For 7jj, the proof is even simpler: <p(A) = A x g\ x er_i is always in T3 
since <j\ x <t_i g 7ij and an equicontinuous point in a Cartesian product 
induce equicontinuous points for each component. Theorem 15 .41 concludes 
for the faithfulness. 

■ 

These encodings allow to transport some properties of general C A concerning 
the top of quasi-orders into order structures induced by each familj|f|. 

Corollary 6.1. Let =4 be or and let $ be a family of CA among: CA 
with two states, T2, T3, captive CA. Then we have the following properties: 

• the set of =4 -universal CA in $ is not co-r.e. 

• for any non-universal A G J, there is a non universal B G $ with A =4 B 
but B/^A. 

Proof. The first property is a direct corollary of theorem 16.41 and 15.31 bv defi- 
nition of faithful encodings. 

For the second property, consider the encoding <f> established in theorem 16.41 
and let A G J be any non-universal CA. By theorem 15.41 there is some non- 
universal CA B such that A =4 B but B/^A. By faithfulness of </>, 4>{B) G 5 is 
not universal and by the definition of encoding it simulates A without being 
simulated by A. I 



8 A stronger result concerning captive CA appears in 37]: ^-universality is undecidablc 
even if we restrict to captive CA with a fixed (but sufficiently large) radius. 
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The families considered above induce structures sharing some properties with 
the general quasi-orders 'near the top'. However, the complete characterisation 
of equivalence classes occupied by some CA of these families is more challenging. 

Open Problem 9. What are the quivalence classes of the simulation quasi- 
orders containing a 2-states CA? a captive CA? a CA from T2? a CA from 
%? 

7. Summary of results 

Figures [3j [5] and [4] hereafter give a summary of results and open problems 
concerning each of the three quasi-orders studied in the paper. 
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Thm. [5"T1 




Fi gure 4: Mixed bulking (quasi- order 
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Open Pb. [7] 
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